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Numerical simulations of dark matter collapse and structure formation show that in addition to a 
large halo surrounding the baryonic component of our galaxy, there also exists a significant number 
of subhalos that extend hundreds of kiloparsecs beyond the edge of the observable Milky Way. We 
find that for dark matter (DM) annihilation models, galactic subhalos can significantly modify the 
spectrum of electrons and positrons as measured at our galactic position. Using data from the recent 
Via Lactea II simulation we include the subhalo contribution of electrons and positrons as boundary 
source terms for simulations of high energy cosmic ray propagation with a modified version of the 
publicly available GALPROP code. Focusing on the DM DM 4e annihilation channel, we show 
that including subhalos leads to a better fit to both the Fermi and PAMELA data. The best fit 
gives a dark matter particle mass of 1.2 TeV, for boost factors of Bmh = 90 in the main halo and 
Bsu = 1950 — 3800 in the subhalos (depending on assumptions about the background), in contrast 
to the 0.85 TeV mass that gives the best fit in the main halo-only scenario. These fits suggest that at 
least a third of the observed electron cosmic rays from DM annihilation could come from subhalos, 
opening up the possibility of a relaxation of recent stringent constraints from inverse Compton 
gamma rays originating from the high-energy leptons. 



1. INTRODUCTION AND SUMMARY 

Recent observations of the spectrum of electrons and 
positrons by the Fermi collaboration [1] and of the 
positron fraction e"'"/(e"'" + e~) by the PAMELA experi- 
ment [2] hint at a possible new source of cosmic ray e+ 
and e~ in the TeV energy region. According to recent 
models [3-9], these excesses could be the signal of dark 
matter (DM) annihilation via a dark sector gauge bo- 
son that allows a Sommerfeld-type enhancement at low 
velocities. Best fits to the electron-positron spectra in- 
dicate that the dark matter candidate x that annihilates 
within the galaxy should have a mass of around M-^m — 1 
TeV, and annihilate into two pairs of light leptons via the 
process DM DM ^ ^ Ae or DM DM (jxj) 4^. 
The particle (j) should furthermore be light enough not 
to decay into pp pairs since excess antiprotons are not 
observed by PAMELA. ^ 

While the visible galaxy spans a diameter of approx- 
imately 40 kpc and a height of 8 kpc, A^-body simula- 
tions [10, 11] predict a roughly spherical structure of dark 
matter subhalos whose peak concentration occurs ~ 70 
kpc from the galactic center (GC) and extends as far as 
several thousand kpc. Relative velocities between parti- 
cles in these regions are one to two orders of magnitude 
smaller than in the Milky Way's main halo, and the rela- 
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^ however see ref. [7] for arguments that this constraint may not be 
necessary, due to astrophysical uncertainties in the background 
model 



five overdensity of such regions make them ideal sources 
of DM annihilation products. This has been explored 
by other authors in the context of gamma ray signals 
originating from subhalos [8, 12-16] and found to be sig- 
nificant. These gamma rays, which originate from final- 
state radiation, are not the main product of this class of 
DM annihilation; rather they are by-products of charged 
particles and neutral pions. 

In this work we consider the possibility that the excess 
leptons observed by PAMELA and Fermi/LAT them- 
selves have a strong component originating in the sub- 
halos. This possibility was previously considered in ref. 
[17], but there it was assumed that one or two nearby 
subhalos would dominate any additional contribution to 
the signal. Here we will show that the best fits to the 
data are found by taking into account the full ensemble 
of substructures. It will be seen that the subhalos that 
individually contribute weakly to the lepton flux are nev- 
ertheless so numerous that their combined effects cannot 
be neglected. 

We used a modified version of the GALPROP cosmic 
ray propagation code, in which leptons from distant sub- 
halos give a new source term at the boundary of the 
diffusion zone. The data of the Via Lactea II simulation 
[10] are taken as our model for the subhalos. We allow 
for independently adjustable boost factors for the main 
halo and subhalos, motivated by the fact that Sommer- 
feld enhancement of the annihilation cross section can be 
much greater in the subhalos due to their lower velocity 
dispersion [3] . If we also allow the background electron 
and positron flux normalizations to be rescaled, as in ref- 
erences [9, 18, 19], we find that the inclusion of subhalos 
gives a much better fit to both cosmic ray data sets, with 
the best-fit DM particle mass of Mom = 1-2 TeV. 

On the other hand, if the e"*" and e~ backgrounds are 
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instead fixed at tlie GALPROP output level, it is known 
that there is a discrepancy between the boost factors 
needed for explaining PAMELA and Fermi, even in the 
standard main halo-only scenario. This discrepancy re- 
mains in the subhalo scenario, where we find that a DM 
mass of Mdm = 2.2 TeV improves the fit to the Fermi 
data, whereas the fit to PAMELA is not improved. 

Our results suggest that the inclusion of leptons from 
DM annihilation in the subhalos surrounding the Galaxy 
should affect not only the amplitude, but also the shape 
of the observed spectrum due to inverse Compton scat- 
tering (ICS) of the leptons with the radiation fields inside 
the observable galaxy. This is potentially important be- 
cause the most recent constraints on this effect [18, 19] 
effectively rule out the DM annihilation interpretation of 
the Fermi excess, and leave only a very reduced corner of 
parameter space consistent with PAMELA. We hope to 
quantitatively address the question of whether substruc- 
ture indeed allows for relaxation of these constraints in 
the near future. 

In section 2 we briefly describe GALPROP and our 
choices of parameters for cosmic ray propagation. Section 
3 details the modifications we made to GALPROP in 
order to include the e~^e~ pairs from DM annihilation in 
the subhalos. The results are presented in section 4, and 
conclusions in section 5. 



2. COSMIC RAY PROPAGATION MODELS 

Inside the diffusive zone of the Galaxy, cosmic ray 
species propagate according to the transport equation 
[20] 

^ = q{r, z,p) + V- {D,,Vi> - KV') (1) 

5 9 „ 5 1, d I dp p ^ , , 

il){x,p,t) is the particle density per unit momentum 
p = \p\, q{x,p) is the source term, D^x is the energy- 
dependent diffusion coefficient, Dpp quantifies reaccel- 
eration via diffusion in momentum space and Vc is the 
convection velocity. In the case of composite species, 
terms accounting for radiative decay and fragmentation 
must furthermore be included. For D^x we use the 
parametrization [21] 

Dxx = Dqxx ^4 Gev) ' ^'^^ 

where E is the particle energy and Dqxx is the diffusion 
coefficient at reference energy E = A GeV.^ Dqxx and 



the exponent 8 are determined by fitting to heavy nuclei 
cosmic ray data. 

There are two widely-used approaches to cosmic ray 
propagation within the galaxy. The first is a semi- 
analytic model in which the baryonic component of the 
galaxy is accelerated in a thin disk at z = from which 
particles diffuse according to a Bessel series expansion 
until z = iicff, beyond which they freely escape. The 
second, fully numerical, approach implemented in the 
publicly available GALPROP [20] package uses a Crank- 
Nicholson scheme to solve eq. (1) within a diffusion zone 
of height LcB and radius i?off. An advantage of the lat- 
ter technique is that it allows the use of realistic maps 
of radiation and gas in the propagation scheme. While 
this is not the focus of our paper, it is relevant to point 
out that differences between models arc responsible for 
differences between fits in recent dark matter annihila- 
tion models. This discrepancy has been known for some 
time; see for example the discussion in ref. [22]. Nev- 
ertheless, we shall henceforth focus exclusively on the 
numerical approach. Our simulations were run using a 
modified version of GALPROP 50. Ip that was graciously 
provided by the authors of ref. [6] and which we further 
modified to handle subhalo sources. 

The strongest available constraints on cosmic ray prop- 
agation models are ratios of secondary-to-primary species 
such as B/C or sub-Fe/Fe. The authors of ref [21] con- 
ducted an exhaustive search of the GALPROP parame- 
ter space for input values that gave best fits to 12 sec- 
ondary/primary cosmic ray experiments. For our simu- 
lation runs we took their best fit parameters: Dqxx = 
6.04 X 10^8 cm2 s-i (0.19 kpc^/Myr), L^g = 5.0 kpc, 
5 = 0.41, with no convection. We used an Alfven speed 
Va = 31 km s~^, which gave a slightly better fit to the 
HE AO B/C data [23]. It should be noted that these 
parameters are quite different from the corresponding 
best fits of the semi-analytic model used, for example 
by Meade et al. [9]. 



3. INCLUDING SUBHALO FLUX IN GALPROP 

The many-body simulation Via Lactea II [10], which 
modeled the evolution and collapse of more than 10^ 
particles over the history of a Milky Way-sized struc- 
ture, resolves over 20,000 dark matter subhalos around 
the galactic host halo. The data characterizing each of 
these subhalos is publicly available [24]. While the visible 
galaxy is only some 40 kpc across, these subhalos extend 
as far out as 4000 kpc from the galactic center. Each 
subhalo is locally much denser than the host halo and 
has its own radial velocity dispersion profile. The anni- 
hilation rate of dark matter, proportional to p^, should 
thus spike within these subsystems when compared to 



^ More precisely, difFusion depends on particle rigidity, the energy 
divided by the charge. We assumed the particles have unit charge 



here. 
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the annihilation rate of diffuse DM particles of the host 
halo. 

For a given subhalo i at a distance £i from the edge of 
the diffusion zone of the galaxy, the flux of e+ or on 
this boundary takes the form 



subhalo i 



2 2 



dr 



(3) 



where Ssh is an average boost factor for the subhalos due 
to Sommerfeld enhancement for example, and pi{r) is the 
mass density profile of the subhalo. The unboosted cross 
section is assumed to be {av) = 3 x 10~^^ cm-^ s~^ in 
accordance with the standard assumption that the DM 
abundance was determined by freeze-out starting from a 
thermal density. In a more exact treatment, the boost 
factor would be velocity dependent [25, 26] and appear 
within the average over DM velocities indicated by the 
brackets in {av). Moreover each subhalo in general has a 
different boost factor since the velocity dispersions that 
determine B^^ depend on the size of the subhalo [13]. For 
this preliminary study, we simply parametrize the effect 
by an average boost factor, where the averaging includes 
the sum over all subhalos as well as the integration over 
velocities. 

The energy spectrum dN/dE of electrons from the DM 
annihilations is taken for simplicity to be a step function 
at the interaction point, dN/dE — M~^^Q{Mj^m — Eq), 
where Eq is the energy immediately following the an- 
nihilation. We are interested in models where the DM 
particles initially annihilate into two hidden sector gauge 
or Higgs bosons, each of which subsequently decays into 
e+e" [3]. The four-body phase space would thus be a 
more exact expression for dN/dE, but the step function 
has the correct qualitative shape and is simpler to imple- 
ment in GALPROP. 

The energy of the electron at the edge of the galaxy is 
reduced from its initial value Eq by scattering with CMB 
photons before reaching the galaxy (starlight, infrared 
radiation and synchrotron radiation are only important 
in the inner galaxy [27]), according to the loss equation 
dE/d£ = -kE^ [9] where k = {Aax/'^ml) ucmb = 6.31 x 
10"'' kpc GeV~ , ctt = ^{aEM^/ '^ecY is the Thomson 
cross-section and ucmb = 0.062 eV/cm'^ is the present 
energy density of the CMB. It is convenient to write the 
solution in the inverted form: Eq — {—kI+1/E{£))^^ 
for substitution into dN/dE. Numerically, we find that 
the losses outside the diffusion zone make a small correc- 
tion, and that the distinction between Eq and E{t) is not 
important here. 

Each subhalo is characterized by a density profile that 
has been fit to the Einasto form 



Ps,i exp 



r 



- 1 



(4) 



with a = 0.17 [15]. The scale radius is found to be 
proportional to the radius Ty^^^^ at which the veloc- 




8 kpc 



20 kpc 

FIG. 1: Geometry of a subhalo shining leptons on the bound- 
ary of diffusion zone of the galaxy. 



ity dispersion is at a maximum, through the relation 
Ts = r„„^^/2.212, while the prefactor scales with the 



maximum velocity Wmax as Ps 



^/(0.897 • 47rr2G). 



To incorporate the contribution (3) to the lepton flux 
from the subhalos in GALPROP, we add delta func- 
tion source terms to q{r,z,p) for the cylindrical surface 
bounding the diffusion zone, as illustrated in fig. 1: 



9disk = 2S{z±h/2)Y, 



d^ 

dE 



<?band 



(5) 



where h and R are respectively the height and radius of 
the cylinder. The factor of 2 corrects for the fact that 
sources in GALPROP have no directionality, whereas the 
fiux impinging on the surface is inward. The sum is over 
the 20,048 resolved subhalos in the Via Lactea II simu- 
lation. In addition, the sources were averaged over the 
aziniuthal angle (p because GALPROP assumes cylindri- 
cal symmetry in its 2D mode (and the 3D mode runs too 
slowly for our purposes). Finally, the distance £i must 
be corrected for subhalos that are close to the diffusion 
zone; rather than the distance to the center of the galaxy, 
it should be the distance to the cylindrical boundary. On 
average, this is a reduction by 17 kpc compared to the 
distance to the galactic center. 

Although most subhalos were located outside of the 
diffusion zone, there are 143 lying inside, whose contri- 
bution required special treatment. Assuming approxi- 
mate isotropy, we took their entire flux of e+ -I- e~ to be 
pumped into the diffusion zone from the boundary rather 
than from their individual positions. Treating them in 
this manner allowed us to group their contribution with 
that of the other subhalos and thus consider a single av- 
erage boost factor for all subhalos. This approximation 
would break down if one subhalo happened to be very 
close to our position in the galaxy, but treating such a 
case would anyway require going beyond the standard 
cylindrical symmetry (2D) mode of GALPROP and us- 
ing the much slower 3D mode. We believe this treatment 
is conservative in the sense that it should only underes- 
timate the contributions of the nearby subhalos. 
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4. NUMERICAL RESULTS 

We compared the observed flux of positrons and elec- 
trons generated by dark matter annihilation within the 
main halo (MH) to a scenario in which both main halo 
and subhalo (MH+SH) DM annihilation occurs. We re- 
stricted our analysis to the DM DM — > 4e channel. This 
channel is simpler to analyze, and is somewhat less con- 
strained by inverse Compton gamma ray constraints than 
the other 4-lepton final-state models, or those with only 
two leptons [9]. 

For the MH only scenario, the annihilation cross- 
section was augmented by a constant boost factor -Bmh, 
representing the effect of Sommerfeld or some other kind 
of enhancement [3, 28-30]. This was varied in order to 
find a best fit to each data set. A similar approach was 
used in the case of MH+SH, where we varied the MH and 
SH boost factors independently. This is justified by the 
expectation that Sommerfeld enhancement should be sig- 
nificantly larger in the subhalos due to their lower veloc- 
ity dispersions [3]. The subhalo boost factor Bsh might 
also have further contributions besides Sommerfeld en- 
hancement, such as the presence of unresolved subhalos 
that we do not take into account [15], as well as substruc- 
ture within the subhalos themselves [13]. 

We minimized the chi squared coefficient 



E 



Ddcl) 



(6) 



where the sum runs over the measured or predicted values 
of ^, which stands for either E^d{ipe+ + 4'e-)/dE in the 
case of the 25 Fermi data points, or V'e+/(V'e+ + V'e-) 
in the case of the PAMELA data, and V'e* is the flux 
of electrons or positrons. When fltting to PAMELA we 
excluded the first 8 of 16 data points, following the usual 
assumption that the dip relative to the background is 
accounted for by solar modulation effects [7]. 



4.1. Freely-varying background 



Freely- varying background 






XPcrmi 


Xpamela 


Xtotal 






MH 

MH-fSH 


0.85 TeV 
1.2 TeV 


15.5 
2.3 


18.7 
14.2 


34.3 
16.5 


90.3 
92.8 


3774 


Fixed GALPROP background 


MH 

MH+SH 


1.0 TeV 
2.2 TeV 


8.2 
2.1 


144 
175 


152.2 
177.1 


110 
146 


1946 



TABLE I: Best fit scenarios. Top: when the background 
positron and electron spectra were allowed to vary by an over- 
all factor; this corresponds to the best overall fit to the data. 
Bottom: using background that was fixed at GALPROP 's 
normalization. In this case we used the best fit to Fermi, 
since the best overall fit gave values of electron + positron flux 
that were ruled out by the Fermi data. MH: main halo DM 
annihilation only. MH+SH: subhalo annihilation included. 
ruDM is the DM mass that gives the best fit and Xi s^re the 
chi squared fits to the respective experiments as described in 
eq. (6). -Bmh and -Bsh are the boost factors necessary for 
MH and SH annihilation cross sections, respectively. Note 
that the addition of a subhalo contribution greatly improves 
the best fit for both Fermi and PAMELA. The required DM 
mass is larger because of the energy loss suffered by electrons 
propagating to us from the galactic edge. 
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Our best fits to the PAMELA and Fermi data were 
obtained by letting the astrophysical background elec- 
trons and positrons be rescaled by overall normalization 
factors, which was also the approach taken in references 
[9, 18, 19]. Adding subhalo contributions significantly 
improved the fits to both the Fermi and the PAMELA 
data. While the best overall fit with only MH electrons 
was for Mdm = 850 GeV (x^tai = 34.3), the MH+SH sce- 
nario gave a best fit at Mr,M — 1-2 TeV (Xtotai ~ 16.5). 
In this case 31% of the DM electron + positron flux at 
the sun's location originated from subhalos. 

A summary of these results is presented in the top 
portion of Table I. The predictions for MH and MH+SH 
scenarios are shown for the total e+ + flux in figures 2 
and 3 and for the positron fraction in figures 4 and 5. The 
value of versus A/dm is shown in fig. 6, marginalizing 



FIG. 2: Fermi data and prediction for e"*" + e~ flux of the 
best main-halo-only fit to Fermi and PAMELA data, with an 
unconstrained background. 



over the background normalizations. For the minimum 
point of the MH+SH model, the background electrons 
had to be reduced to 97% of their predicted values, while 
background positrons were rescaled to 137% of the GAL- 
PROP output. 

The optimal boost factors of 90 for the main halo and 
3800 for the subhalos are quite reasonable from the point 
of view of DM models that give Sommerfeld-enhanced 
cross sections [3]. We leave for future work the issue of 
detailed particle physics model building to match these 
and other features of the best-fitting models. 
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FIG. 3: Same as fig. 2, but now including subhalo contribu- 
tions to the lepton fluxes. 
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FIG. 4: PAMELA data and predicted positron fraction of tlie 
best main-halo-only fit to Fermi and PAMELA data, with an 
unconstrained background. 
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FIG. 5: Same as fig. 4, but now including subhalo contribu- 
tions to the lepton fiuxes. 



Subhalo and main lialo anniliilation 

IVIain lialo only 




DIVI mass (TeV) 

FIG. 6; Combined for the Fermi and PAMELA data as 
a function of the dark matter mass, for the unconstrained 
background. Dashed (blue) line: main halo DM annihilation 
only. Solid (red) line: subhalo and main halo contributions 
combined. 



Subhalo and main halo annihilation 

Main halo only 




DM mass (TeV) 

FIG. 7: versus Mom for the Fermi e"*" -I- e~ data using 
the GALPROP constrained background. Dashed (blue) line: 
main halo only. Solid (red) line: subhalos plus main halo. 



4.2. Constrained background 

We performed a second analysis by taking the elec- 
tron and positron backgrounds to be those predicted by 
GALPROP. In this case, although there is no good simul- 
taneous fit to the combined PAMELA and Fermi data, 
we nevertheless find that SH contributions improve the 
fit. In rough agreement with ref. [6], we find that the 
PAMELA data require a boost factor several times higher 
than that needed to fit the Fermi data. 

The plots of versus M^m, for both the MH-only and 
MH+SH models, are shown respectively for the Fermi 
and PAMELA data in figures 7 and 8. It is striking that 
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— Subhalo and main halo annihilation 
- Main halo only 




DM mass (TeV) 

FIG. 8: 'X versus Mdm for the Pamela positron fraction data, 
using the GALPROP constrained background. Dashed (blue) 
line: main halo only. Solid (red) line: subhalos plus main halo. 



the best-fit DM mass becomes significantly larger and 
less constrained in the fit to the Fermi data including 
subhalos, fig. 8. The increase in the required DM mass is 
due to the energy lost by electrons and positrons during 
propagation from the edge of the diffusion zone to our 
position. The best fit to the Fermi data has AfoM = 2.2 
TeV (x^ = 2.05) with SH+MH, compared with A/^m = 1 
TeV (x^ = 8.15) in the MH only case. The required boost 
factors for these fits are Smh = 146 and i?sH = 1946 for 
SH+MH, in contrast with Smh = 110 for MH only. The 
best fit cases for the e^+e~ spectrum are shown in Figure 
9 and the results are summarized in the bottom part of 
Table I. The corresponding positron fraction in each of 
these scenarios is shown in Figure 10. 

The fit to PAMELA is also improved by the addition 
of SH positrons, but only at low DM mass, < 500 

GeV. However the best fit parameters for the PAMELA 
data by themselves lead to a prediction of the e+ + e~ 
E^dN/dE spectrum that exceeds the Fermi data by more 
than 3cr, resulting in a = 460 fit to Fermi. The badness 
of this fit is evident in fig. 11. 



4.3. Relative contributions of subhalos 

It is interesting to quantify how much of the signal can 
be contributed by the subhalos relative to that coming 
from the main halo. We show the fraction of e++e~ pairs 
due to the subhalos, as a function of the DM mass, in fig. 
12. For the best-fit values of the mass, this fraction is 
around 30%, but for larger values of Mom (yet still giving 
reasonable fits) it rises to 60% or more. This may be 
helpful for weakening the constraints on the model from 
production of gamma rays by inverse Compton scattering 
[18, 19]. We hope to investigate this issue in the near 



10 - 



> 

a> 
O 




Total flux 

Astrophysical 

Main Halo DM contribution 

Fermi data 



10 



10 



Energy (GeV) 
(a) 




Main Halo Contribution 

Astrophysical Background 

Fermi Data 



10 



10 



Energy (GeV) 

(b) 

FIG. 9: Best fits to Fermi data, a) Main halo only, Mdm = 1 
TeV, with a boost factor Bmu = 110. b) Subhalos plus main 
halo. A/dm = 2.2 TeV and Bmh = 146, Bsh = 4825. 



future. 

Another relevant issue is the hierarchy of contributions 
of subhalos relative to each other. One would like to 
know whether it was really necessary to add the contri- 
butions of all 20,000 subhalos, or if perhaps only the few 
closest ones dominate. Fig. 13 shows the distribution of 
subhalos contributing a given flux $ (normalized to the 
contribution of the subhalo that gives the largest value 
'I'lnax), weighted by the flux, and also the integral of this 
quantity. From the integral, we see that 50% of the total 
signal comes from subhalos whose individual intensities 
are less than 5% of the strongest one. Thus to get a quan- 
titatively accurate estimate, it is necessary to include the 
very numerous subhalos whose intensity is low. This also 
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FIG. 10: Positron fraction for four of the best fit scenarios 
with constrained backgrounds. Top (sohd) hues correspond 
to fits to PAMELA data only. Uppermost (magenta): MH 
only; tudm = 250 GeV, Smh = 225. Lower (red): MH+SH; 
moM = 150 GeV, Smh ~ 9.3, Ssh = 509. Bottom (dashed) 
lines correspond to the best fits of these scenarios to the Fermi 
data. Upper dashed (blue): MH only; moM ~ 1 TeV, Smh 
= 110. Lower dot-dashed (black): MH+SH; moM = 2.2 TeV, 
Smh ~ 146, Ssh = 4825. Although the former set provide 
a better x^, they predict a total e"*" + e~ flux that conflicts 
with the Fermi data by at least 3cr (see flgure 11). 
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FIG. 11: e+ + e" curve for best flt to PAMELA data only. 
niDM = 150 GeV, Smh = 9.3, Ssh = 509. Although the 
high energy tail could be compensated by other sources (e.g. 
pulsars), the fact that the model exceeds the data points at 
low energy leads us to disfavor this model. Note that the 
subhalo contribution is too small to be seen in this figure. 



FIG. 12: Proportion of the total fiux of e + e~ originating 
from subhalos as opposed to the main halo, as observed 8.5 
kpc from the galactic center (the position of the solar system) 
in order to obtain a best fit to the Fermi data. Each point 
represents an individual simulation. 



5. CONCLUSIONS 

We have shown that the inclusion of electrons and 
positrons from the galactic subhalos can significantly al- 
ter the predictions from annihilating dark matter models. 
Using the Via Lactea II simulated data of the subhalo 
distribution around a Milky Way-like galaxy, we found 
that the contributions from substructure can give im- 
proved fits to the PAMELA and Fermi excess lepton data, 
and increase the value of the expected mass of the dark 
matter particle. A strong Sommerfeld boost coming from 
the low velocity dispersions of the subhalos, as well as 



Ir 




max 



suggests that our computation is an underestimate, since 
we do not count the subhalos that are not resolved by the 
Via Lactea II simulation. 



FIG. 13: Vertical bars (black): $diV/d$, the distribution of 
subhalos contributing a given fiux at the edge of the difusion 
zone, weighted by the flux. Continuous curve (red): integral 
of $ dN/A^. 
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the uncounted contributions of subhalos unresolved by 
the Via Lactea II simulations, are possible sources for 
the boost factor necessary to obtain best our best fits to 
the data. According to these fits a third or more of the 
electron cosmic rays from DM annihilation could come 
from subhalos outside of the visible Milky Way. 

The next step for future work will be to see whether 
the reduction of the flux from the main halo can weaken 
Fermi constraints on annihilating DM models due to 
the inverse Compton gamma rays produced by the high- 
energy leptons [18, 19]. These constraints are sufliciently 
strong to rule out the DM interpretation of the Fermi lep- 
ton excess, under the usual assumption that all the e+e~ 
pairs are produced in the main halo. The constraints are 
strongest from data near the galactic center. By shifting 
the production away from the center to the subhalos, the 
constraints should be weakened, but whether the effect 
is large enough to reinstate the DM interpretation of the 
Fermi lepton observations is a quantitative question. In 
addition, one should satisfy other protohalo constraints 



[31], extragalactic gamma background [32] and last scat- 
tering surface CMB constraints [33]. 

If it is possible for the scenario to pass these tests, 
it will be interesting to check whether specific particle 
physics models are able to give the average boost factors 
that we have treated as free parameters in this prelimi- 
nary study. 

As we were completing this work, ref. [34] appeared, 
which presents an analytical method for taking into ac- 
count the efl^ect of substructure on dark matter annihila- 
tion. 
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